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ABSTRACT 


A detailed a^aalysis of temperature variation withm 
the sodium-sodium intermediate heat exchanger of the fast 
Breeder test reactor at Kalpakkam has been performed This 
analysis avoids the use of empirically determined heat 
transfer coefficients since it may cause large errors in 
the case of liquid metals The analysis involves a two- 
region Sturm-liouville problem with both negative and positive 
eigenvalues The results of the present analysis predict 
the efficiency of the heat exchanger to be six percent less 
than that predicted by the approximate method using Musselt 
numbers for constant heat flux* It is shown that the 
approximate method can lead to substantial errors if the 
ratio of heat capacity rates of the two streams is around 
2 (or 0 5) and the heat exchanger is not long 



CHAPTER 1 


INTRO rUCTIOR 


[11] INTRODUCTION 

Nuclear power competes commercially at present with 
any other conventional source of power like themal power or 
hydel power There are more than hundred nuclear power plants 
in commercial operation, at present, all over the world 

Nuclear power is produced by controlled fission of 
Nuclear fuel like U'ranium~255 or Pluto mum- 239 Each fission 

releases energy which is responsible ultimately for generating 
power This chain reaction is carried in a reactor core The 
core IS normally composed of fuel, moderator, coolant, control 
devices and structural material Coolant which carries heat 
from hot fuel elements may also act as a moderator in the 
core If moderator is removed from the core and fission 
chain reaction is maintained without slowing down the 
neutrons, then the reactor is called fast reactor Its 
counterpart with moderator is known as thermal reactor 

Transfer of heat to the coolant inside core, is mainly 
by conduction and convection This coolant is pumped to an 
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intermediate heat exchanger, (IHX) A secondary fluid picks 
up heat in IHX and goes to steam generator (SG) Steam produced 
in SG runs steam turbine which is coupled to electric genera- 
tor Therefore the major difference in conventional and 
nuclear power plant lies in the mode of generation of primary 
heat 

[1 2] fast EE actors 

In fast reactors the neutrons are not deliberately 
slowed down by a moderator but only by inelastic collisions 
of neutrons with fuel and other reactor material nuclei, Pig (1 1 
Fast reactors have very high power densities. A typical fast 
reactor has power density about 40 times that of heavy water 
moderated natural uranium thermal reactor 

Liquid metals are the only possible coolants for fast 
reactors, because of high power generation per unit volume, 
within the core. Liquid metals have excellent heat transfer 
characteristics which results in low hot spot factors, i e , 
less local temperature variations Use of liquid metals 
give high thermal efficiency* 

There is very hi^ neutron flux inside reactor core 
This results in hi^ level of induced radioactivity in the 
reactor materials and liquid rdetal coolant Coolant is 
circulated in the heat exchanger and therefore there is q, 
chance of leakage of radioaotive coolant to the other ei^e of 
the heat exchanger 
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If on the other side, the working fluid is steam then 
the leakage of radioactive coolant to fche steam side may 
expose the wo iking personnel (dealing with the turbine 
or condenser) to harmful radiation 

Radio active sodium has a long half life of 15 hrs 
and it emits high energy gamma rays Hence leakage of 
sodium to steam cycle is dangerous This is avoided by 
employing an intermediate heat exchanger (IHX), in which 
liquid sodium flows on both the sides This reduces the 
possibility of radioactive contamination of steam because 
now the steam is generated in steam generator, where the 
secondary sodium (not active) is flowing on the other side 

[15] PAST BREEDER REACTORS (PBR) 

Interest in fast reactors stem from the fact that 
they offer a singular possibility of converting fertile 
material to fissile material Past reactors can have ratio 
of conversion of fissile material generated (from fertile 
material) to fissile material consumed, equal or greater 
than one Such a reactor is termed as fast breeder reactor 
(PBR) 

PBR thus offers the possibility of producing power and 
nuclear fuel simultaneously This can be achieved economi- 
cally only when many of the technological hurdles are 


overcome 
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[ 14 ] LIQUID MDms 

Past reactors do not have moderators and therefore it 
must be ensured that coolant does not act as moderator Hence 
the coolant nucleus must have small cross-section for inte3>- 
action with neutrons This coolant has to satisfy stringent 
requirements as regards its thermal properties As mentioned 
earlier the rate of heat generation inside fast reactor core 
IS very high Therefore coolant must have high thermal 
conductivity, low vapour pressure, high specific heat and 
high volumetric heat capacity These properties will enable 
it to transport heat from core rapidly with less pumping 
power and smaller area of contact For reactors with high 
thermal fluxes and operating at high temperatures, liquid 
metals have special interest as coolants Liquid metals 
have excellent thermal properties and they fulfil the 
above mentioned requirements quite well 

Amongst various liquid metals sodium is the most 
suitable coolant for a reactor operating at high temperatures, 
particularly for fast reactors where moderating (low mass 
number) materials cannot be used In fact sodium and sodium 
potassium alloys have been successfully used as coolants in 
both fast and thermal reactors^ They are also used as 
secondary coolants in intermediate heat exchangers, in 
nuclear reactor power plants 
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[15] LlQUIIi METAIi HLAT EXOHANGERS 

Research in liq^aid metal heat transfer area is going 
on for last three decades Most of it was devoted to develop- 
ing correlations for use in engineering design Accurate 

and reliable prediction w-exe not considered important till 
recently The main reason for this was relatively poor 
thermal properties of heat exchanger wall material This 
resulted in the temperature differences across the wall 
thickness an order of magnitude higher than temperature 
differences associated with convection to liq,uid metals 
Therefore, fairly large inaccuracies (say + 30 percent) of 
liquid metal heat transfer coefficients did not greatly 
affect the overall thermal design But for liquid metal 
heat exchangers used in fast breeder reactors a very accurate 
thermal design is required on account of safety considerations 

Liquid metals have certain peculiar properties In 
liquid metal heat exchangers the use of conventional heat 
transfer coefficient concept is not valid, hence the existing 
methods for nonmetallic fluids are not directly applicable 
to liquid metals 

Liquid metals have small viscosities Therefore 
their application as coolant or heat exchanger fluid, e g,, 
in nuclear reactors, almost always involves large Reynold 
no and consequently turbulent flow, Thermal conductivities 



6 


of liquid metals are very high, hence molecular conduction 
remains significant throughout the cross-section of flow 
channel or duct geometry Therefore temperature gradients 
are not localized near the wall as it is the case with 
nonmetallic fluids This is a characteristic difference 
between turbulent convection in liquid metals and nonmetals 

Prandtl numbers for liquid metals are very small 
Therefore, Peclet numbers associated with even turbulent 
flow may be as low as 50, with nonmetallic fluids it is 
rarely so Low Peclet number associated with high thermal 
conductivity results in axial heat conduction Millsaps 
and Pohlhausen (1956) and Singh (1957) have shown that upto 
a Peclet number of 100 axial conduction is significant But 
error due to neglect of axial conduction even at a low 
Peclet number may not be large Schneider (1956) determined 
that for Peclet number of 10 and for length to diameter 
ratio of unity, the neglect of axial conduction gave error 
of only 4 percent 

Liquid metal Husselt numbers are very sensitive to 
boundary conditions Sleicher and Tribus (1957) have 
computed the ratio of Fusselt numbers for uniform heat 
flux to that for uniform wall temperature for different 
fluids as a function of Eeynolds number This is shown 


NUSSEU NUMBER RATIO 


FIG 1 2 



Effect of Dounc conditions on fully 
developec Nuss , numbers; ratio of 
values fo unite ornr 

ter 
(V^ 


c 
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in Pig (1 2) Prom this figure -vre conclude that non 
metallic fluids are insensitive to the type of boundaiy 
conditions specified Metallic fluids on the other hand 
are very sensitive to this and ^variation in Nusselt 
numbers for two different ca‘='es may be as high as 40 percent 

In liquid metal heat exchangers the heat flux is 
rarely uniform on the other hand boundary conditions which 
specify the wall temperatures are not often accurate 
description of actual situations of practical interest in 
liquid metal heat transfer In certain simple condensers 
or boilers^ wall temperature boundary conditions can serve 
as reasonable approximation^, but with liquid metals this 
approximation may be quite inaccurate* Por example, because 
of excellent heat transfer properties of liquid metals , 
conduction through the condenser or boiler walls will be 
much more important than with nonmetallic fluids 

With heat transfer coefficients sensitive to boundary 
conditions and large number of boundaiy conditions that are 
of practical interest, most of the convenience of traditional 
heat transfer coefficient concept is lost Hence one has to 
start from the statement of basic energy equation for the 
given configuration coupled with the careful statement of 
boundary conditions. Assumption of uniform heat flux of 
uniform wall temperature may lead to significant inaccuracies 
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The thermal modelling of IHl of -PETR at Kalpakham 
is considered in this thesis. The temperature profile 
within the heat exchanger is obtained by solving the basic 
differential equations of heat transfer The use of the 
concept of heat transfer coefficient is scrupulously 
avoided This is because it has been shown by btein ( 1966 ) 
and others that the use of heat transfer coefficient 
concept can lead to great errors in the prediction of 1:6111“ 
perature profiles, in liquid metal heat exchangers 
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MOIELLING OP INTBEMEDIASP HEAT EXCHAiaSK AT 

ICAIPAKZAM 


[2,1] IHTBDDUOTIOH 

A fast breeder reactor (PER) power plant is coming up 
at Kalpakkam, near Madras This plant is expected to 
produce power, commercially, from early ninties Its design 
IS based on the Prench: PBR called RAPSODIE a detailed 
description of this reactor system is given by Yevick (1966) 
This PBR can be used to breed Plutonium - 239, a f^-ssile 
material, using Uranium - 238 as its parent fertile material 

A prototype of the commercial reactor is likely to come 
into -- operation in early eigthies Its designed capacity 
IS 45 MW th A detailed description of its intermediate 
heat exchanger is given here 

The prototype is called fast breeder test reactor 
(PBTR) PBTJR IS a sodium cooled reactor The heat produced 
at the reactor core is removed by circulating sodium through 
the core The primary heated sodium then transfers heat to a 
secondary sodium at This secondary sodium then heats 

up the water in steam generator, for running steam turbine to 
produce electric powers 



SECONDARV bODILM iNl 
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There are two IHX in FBTR and they share the load 
eq^ually viz 22 5 MW th Their normal operating condition 
IS given in Table 2 1 


Table 2 1 Normal Operating Condition (For single IHX) 


Power exchanged 
Primary fluid 
Primary fluid flow rate 
Primary inlet temp 
Primary outlet temp 
Secondary fluid 
Secondary inlet temp 
Secondary outlet temp 
Pressure (inlet to unit) 
Pressure (outlet to unit) 
Secondary flow rate 


22 5 M¥ th 
Sodium (active) 

125 kg/sec 
5 20°G 
380°0 

Sodium (not active) 

2d3°0 

514^0 

5 03 bars abs* Primary 

4 81 bars abs fluid 

*74 kg/sec 


[ 2 2] INTEhMEPIATB HEAT EXCHANGER 


There are two IHX They consist of a fixed shell, a 
removable portion having the tube bundle and biological 
shield and other provisions for IHX instrumentation and 
control, Fig (2.1) 
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a) Primary Side 

The primary sodium system components have been 
designed in such a way that the main equipments, e g , heat 
exchanger pumps could be removed without disturbing the 
primary circuit Hence the components have a fixed tank 
connected with the primary circuit in which free sodium 
level IS maintained with Argon gas cover (^ee Fig (21)) 

The removable portion is introduced into this from the top 
The fixed shell of the IHX is a stainless steel vessel The 
vessel has a primary sodium inlet nozzle on its vertical 
portion and outlet nozzle on the dished nozzle enclosure 
at the bottom The vessel has double envelop till the 
height upto which sodium is likely to be present Hot 
nitrogen could be circulated throu^ the nozzle arrangement 
around annulus through the shell and its envelop, to preheat 
equipment before charging with sodium 

Electrical resistance type heaters are mounted on 
double envelope of this vessel for heating up this equipment 
during the initial testing of the secondary loop 

b) Secondary Side 

There are 888 stainless steel tubes of 12/14 nim in the 
tube bundle welded to tube sheets at each end Tubes are 
arranged over 12 pitch circles around a central pipe of 
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diameter 224/254 mm ( see Table 2 2) This central pipe is 
also welded to the tube sheets at either end The secondary 
sodium inlet pipe of diameter 200/208 mm passes inside this 
pipe and is welded only to bottom sheet The annulus between 
two pipes IS kept under argon atmosphere The differential 
expansion between the central pipe and the tubes and also 
among tubes is accommodated by giving a wave shape (expansion 
bends) to the tubes at the top 

There is a vertical baffle around the tube bundle 
There are inlet and outlet windows for primary sodium flow, 
in this vertical shell The shell is free to move vertically 
during thermal expansion By means of a ring provided around 
this shell, it rests over another ring of fixed shell, thereby 
avoiding by pass of the heat exchanger by primary sodium flow 
There are also a few vertical and circumferential baffles on 
the outer side of the vertical shell arranged around the top 
inlet windows to achieve uniform flow distributions around 
windows and also flow without any gas entrainment problem 


Table 2 2 Arrangement of tubes 


Boy 


of tubes 

Rq-W 

No... ..of tuh.es 

1 


41 

7 

77 

2 


47 

8 

83 

3 


53 

9 

89 

4 


59 

10 

95 

5 


65 

11 

101 

6 


71 

12 

107 

Total No of Tubes 888 , Pitch, 

pitch, circumferential 20 mm 

radial 

19 mm 
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There are temporary block devices to prevent incidental move- 
ment of this vertical shell during handling of this eq.uipment 

The entry header for the secondary tube bundle is 
formed by a dishe^d" end welded to the top tube sheet Secondary 
sodium enters at the top, flows all the way down through its 
inlet pipe It collects at the bottom dished end before 
passing through tubes It leaves IHX at the top 

The primary sodium enters throu^ the top of the nozzle 
of the fixed shell It gets distributed around in annulus 
It rises above the nngs described earlier and enters the 
shell through windows at the top It flows downwards from 
there around the vertical tubes which carry secondary spdium 
It finally flows out through the outlet steam nozzle Tb^e top 
dished end collectors and outlet pipe also have electrical 
resistance heater provided for preheating 

c ) Biological Shield Plug 

The plug and the integral part of the removable por'^ion 
of IHX form the closure of this equipment at the top It rests 
over the top flange of this fixed shell and the sealing towards 
the outside atmosphere is assured by means of two '0* rings 
The parts of the heat exchangers which come in contact with 
sodium are made of Austenitic stainless steel AISI 316 
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[2 5] MODJiIjLIN(x 

The IHX IS a shell and tube type heat exchanger It 
IS a single pass counter-flow heat exchanger Secondary 
sodium flows all the way down to the dished end and then 
rises to top exit There is practically no heat transfer 
from or to the secondary fluid during its downward flow 
because of argon shield around it Primary fluid enters 
the annulus around removable shell at a height such that it 
IS close to top entrance windows Therefore this heat 
exchanger can be approximately considered to he a single pass 
counter flow, ahell and tube type heat exchanger Hot fluid 
(primary active sodium) flowing outside tubes and secondary 
fluid (sodium) flowing inside the tubes 

The geometrical configuration of the heat exchanger 
tubes IS such that the problem is one of heat transfer across 
tube bundles for in line flow of fluid This problem has 
been studied by many investigators, employing different 
methods 

a) Equivalent annulus approach was adopted by Dwyer and 
Tu (i 960 ), Priedland and Bonilla (196l) end others 

b) Graphical and lumped parameter methods were used by 
Deissler and Taylor (1957), Dwyer (1966) and others. 


hlG 2 2 (a) 

All dimensions are in mm 



ti 

Q 


Characteristic flow region 

I 

♦ 



FIG 2 2 (b) 



PRIMARY SODIUM 

12 X 14 


Cl ^■■‘eristic ' n converted 

to i/alent 
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c) Finite difference method was employed "by Bender and 
Magee (1969) and Oberjohn (1970) 

[2 4J iiQCJIVALhffl} ALWULU3 APPROACH 

In a shell and tube type heat exchanger, the tube 
centres usually make a finite number of equilateral arrays, 
across the heat exchanger cross-section They are individually 
called unit cell or lattice The cell may be a regular 
polygon eg it is rectangular in IHX considered here, 

Fig ( 2 2a) Bach unit cell represents typically the thermal 
behaviour of entire heat exchanger 

Once unit cell is recognized, characteristic flow 
region is determined around each tube. Pig (2 2a) It is 
also of same gaometry as its parent cell geometry The 
boundaries of the characteristic flow region are considered 
to be adiabatic If this characteristic flow region is 
circular, the mathematical analysis becomes very much sim- 
plified This is unfortunately not the case with almost all 
shell and tube type heat exchangers In equivalent annulus 
approach this characteristic flow region is replaced by an 
imaginary circular characteristic flow region around the tube, 
Pig (2 2b) Outer boundary of the annulus is considered 
insulated The annular flow area is taken equal to the flow 
area of the characteristic flow region in Pig ( 2 2a) 
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Nijsing and Eifler (1973) compared results of 
different approaches, mentioned above, along with the 
approximate analytical method they used They have shown 
that when in equivalent annulus approach the ratio of outer 
radius of annulus to outer radius of tube i e , p ^ 1 15 , 
the results obtained using equivalent annulus approach do 
not differ significantly from that obtained by any other 
method 

The Kalpahkam IHX has p =1 57, therefore we adopt 
equivalent anniilus approach The mathematical simplificatq.O’' 
offered using this method outweigh the insignificant errors 
contained in the final results 




fe double pipe countercurrent flow heat exchanger 


CHiPTER 5 


mathematioal analysis 

[3 1] INTRODUOTIOJ^ 

Employing equivalent annulus appieach the shell and 
tube type IHX at Kalpakkam is converted to an equivalent 
IHX where each tube is surrounded by an imaginary annulus 
Outer boundary of the annulus is adiabatic Hence the 
problem to be solved is c ounter -cur rent , double pipe heat 
exchanger where on both sides liquid metal (sodium) is flow- 
ing A line diagram, Eig (3 1)» represents typxcal annulus 
where in regions 1 and 2 secondary and primary sodium are 
flowing Brown (1968) presented mathematical analysis of 
such a system He basically extended the work of Stein (1S66)* 
Stein discussed the problem in generality considering symme- 
trical and asymmetrical ducts, co-current and counter-current 
flows, with constant and varying temperature and heat flux 
conditions We will use his procedure to formulate the 
problem Final results are obtained using Brown's modifica- 
tions which simplify the numerical calculations 
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f3 2] ASSUMPTIONS 

TNe general equation tor energy balance for double 
pipe beat exchanger precludes any possibility of solving 
it by analytical means Pollowing assumptions are made to 
simplify the equation 

i) The fluids enter the heat exchanger at uniform 

temperature 

II ) Physical properties are independent of temperature 

III ) Frictional heating is negligible 

iv) Axial heat conduction within the heat exchanger 

walls IS negligible 

v) Axial heat conduction within the heat exchanger 
fluids IS negligible 

VI ) Velocity distribution of the fluids is independent 
of axial position (i e , fully developed) and is known 
function of radial variable 

vii) The heatEexchanger is in steady state operation 

The first four assumptions are always almost attainable 
in a physical system Fifth assumption for liquid metals 
IS valid if Peclet number is greater than 50 The sixth 
assumption is realized in liquid metals because of their 
small viscosities They attain fully developed velocity 
profiles very early compared to thermal fully developed 
profiles For large Reynolds number (He No >• 10^) the 
velocity profile is nearly flat even close to entrance region 
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[5 5] DlitoiNblOM POEMUlATIOir 

The general energy equations for tute and annulus, 
with the stated assumptions can be written as follows 


Tube 




[(1 + 



3T^( r^,z) 




1 




3T^( r^, z) 
3 Z 


Annulus 



[(1 



dTJr^,z) 

\ p .1 J 


a. 


-U2(r2) ) 


aT2(r2,z) 

ai 


(3 1) 


(3 2) 


where , 

= Eddy diffusivity of heat in region i 
= Thermal diffusivity of fluid in region, i 
1 = 1,2 

The above equations are valid both for laminar and turbulent 
flow 

Boundary conditions for these equations are as follows 


H 

H 

O 

~ Tiq ( constant ) 

(S 5) 

T2(r2, 1) 

~ *^2L (constant) 

(3 4) 
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Tube centre line, 


outer onou-iits 


ST^(0, z) 




= 0 


3Ti C I 


= 0 


continuity of heat flux across inner wall 


(3.5) 


^12 ^1 




ar-, 


^21 ^2 


9T2(r2i, z) 


sr. 


(3 6) 


conduction of heat across inner wall 


^12 ^1 




w 


3r 


1 In (3:'2 i/^12^ 


(3-7) 


where , 

= thermal condu-ctivity of fluid in region i 

1 s 1, 2 

= thermal conduptivity of tube metal 
[5 4] DliVlENSIONLBSS BO EMULATION 


We introduce following no n-dimensional variables, 


"■l 

X'l ^ 

=s r: — where 0 < rT 

ri2 - 1 

< 1 

(3 

8) 


^2 ^21 

where 0 £ 


(3 

9) 

^2 

^22 ” ^21 

x* 2 

A 

rtr 

4z 

• where 0 ^ 

A 

z 1 L 

(3 

10) 

z 

^®1^ 2^12^ 
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**'/'*' V ^ *• 

1 (t ,z) = i Where 0 < T <1 

rp _m “ 1 ~ 

1 = 1 , 2 (3 11 ) 

and Pe^ = Peclet number of fluid flowing in region 1 

Substituting these parameters in Eq (3«l) and (3 2), the 
governing energy equations become , 


Tube 


~ [f^(r^) r^ 


3r. 


1 


S^TlCri, z) 


] = g3,(?i) 


Annulus 
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(3 12 ) 


T^+a 8 


X- [f2(r2)(r2+a) 


8 T 2 (r 2 , z) 

A 

aro 


= - «2 (- 2 ) - 


IT 


(3 13) 


and boundary conditions become 
Entrance 

5l (r^, 0) = 0 (3 14) 

T 2 L) = 1 (3 15) 

Interior 

z) 

a rj_ 


0 


(3 16 ) 
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where, 


f 


K 


3\(1, z) 




3^2(0, z) 


ar. 


K, 


aT^(l,z) 


w 


3r^ 


z) 


= T.(0, z) - T, (l,z) 


0 


3r. 


(3 16) 


(3 18) 


(3 19 ) 




u 

1 

R 


lav 


H 




K 


= 1 + 


-HI 

a. 


u^(r^) 


U 


lav 


(3 20 ) 
(3 21 ) 


= Dimensional velocity of fluid in region 1 
= Average velocity in region 1 , 


= 1, 2 


^ 21^^22 

R/(l-R) 

*^2^2 


(3 22 ) 
(3 23) 

(3 24) 


cs specific heat of fluid flowing in region i 


= Mass flow rate in region 1 


k 


1 Uzal 

R 


(3 25) 
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In 


(r^i/r^g) 


(3 26 ) 


,2 HKR 

“ l+R 


(3 27) 


[3 5] SOLUTION OT GOVERNING EQUATIONS 

The Eq^ (3 12) and (3 13) are solved using separation 
of variables method This leads to the Sturm-Liouville 
problem for two regions involving discontinuity at interface 
Solution of such problem involves positive and negative 
eigenvalues and their corresponding eigenfunctions A 
detailed mathematical analysis is given in Appendix A, for 
the solution of the Bq (3 12) and (3 13) with the stated 
boundary conditions 


[36] temperature distribution 

The general solutions for temperature distribution, are 
as follows 


Tube 


A A A 


O A ^ 
(1“2) 


= ^ I Ui) e 

n=l 

2 A 

+ e"“ '] (5.28) 


Annulus 


T2(r2,z) 


G + 
0 


+ C 


n 


nil 

^2,n (R' ® 



1 2 ^ 
A Z 

n 


e 


(3 29) 

i 


5 


] 
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where , 







n-th positive eigenvalue 

Absolute value of the n-th negative 
eigenvalue, 

Normalized eigenfunction in region i 
associated with n-th positive eigenvalue 
Normalized eigenfunction in region i 
associated with n-th negative eigenvalue 


G 


n 


Expansion coefficient associated with, 
n-th positive eigenvalue 




Expansion coefficient associated with n-th 
negative eigenvalue 


[3,7] El GENVAEUJji EQUATIONS 

2 2 

The g eigenvalues and are the squares of the 

positive roots of the equations, 


2^ = 0 


2^ (Pn) = 0 


(3 30 ) 
(3 51 ) 


where , 




da(0, x^) dF(l, X„) dG(0, x^) 


n 


n' 


dr. 


A 

dTn 


dr. 


-a-Fd. x^) 


— &(0, X^) 


A 

dr. 


- K 


(3 52 ) 
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= 


, - K 

n w 


dGl:o, 6a(0, 


n' 


dr. 


A 

dr. 


+ 


rj^. 


dr. 


dF’‘(l, p ) . 

P*(l, Pjj) - K <0. P„) (3 53) 


dr. 




1 


The functions S', P*, G-*and G- are defined by follo-wing 
initial value equations for auxiliary functions 

[3 8] AUXILIARY PUNCTIOKS 


, ^ ^ dF(r. , X^) n ^ A 


dr. 


dr. 


(3 34) 


where , 


dF (0, x^) 

^ = 0 , F (0, X^) = 1 


dY. 


dr^^ 

(5 35) 

= 0 , F*(0, p^) = 1 


dr 


dF^O, p^) 


dr. 


[f2 (r2) (rg+c) 


dr. 


dr. 




( 1*2 + a) GCrg, X^) 


0 


(3 56) 
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dr. 


= 0 , 0(1, x^) = 1 


d . * ^ dG*(r-, 6) . . 


dr. 


dr. 




(^P + a) a" (r^, pj,) = 0 


(3 37) 


dG*(l, p^) ^ 

r~^ = 0 , G* (1, p ) = 1 

dr. ^ 


[3,9] BIGBNPUNGTIONS 

The normalized eigenfunctions are given by the following 
equations, 

d&(0, X„) 

"“r.n = C ^ ^ (^1- V]/'rluJ (3 37) 


(^l) 


acfCo, B),. 

[ r— ^ S' (ii- Pn)]/YlUnl <5 58) 


dr. 


B. (r.) 
2,n ' 2^^ 


« (1. J 

[K ;r-^ e <^2, A^)]/ir|Hj (3 39) 


A 

dr. 


^2,n ^^2^ ~ 


djP*(l, Pn^ 


dr. 




where and are defined by following equations for 
normalizing factors 
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[3 10] mmujizim PACTOfito 




n 


[- -:r^] ([BCr^, X„)]^> 


dr. 


P dF(l, X ) p ^ P 

. hk2 [ ^ ]2 ([GKJj, X^)]2) (3 41) 


dr. 


and 


< - C 


dG*(0, 


dr. 


r Bt {[P’‘(r^, 


P dF*(l, 6 ) ^ 5 

HE^ [ _2^j2 B2{[G'‘(r , P„)]^>(3 42) 

dr, ^ n. 


'2 '"2^ ^n' 

where n =; 1,2,3» and j- are defined by Sq. (/ 35) 

and (a 36) as bulk average temperatures over cross-section, 

[3 11] IXPAFSION OOSFFIOIhJNTS 

The expansion coefficients Aj^ and 0^ are defined by 
the following set of infinite linear algebraic equations 


IcL. ®l,k^^l^^ 


2^ 


^ \ 


2 


- % 'Bl,n <^l'l %,k ) °k3 




(3 43) 
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'\.n (r\)> (1- 

®1 ‘®l,n %,k (?i)> (1-e 


® 

-X 2£ 

" )o. ] 


= - ®1 “l,n (%» 

with n = 1,2,3, 

and for the zeroth order coefficient 




A. 


2^ 

“l,k (d)> (l-e"^ ) Cj^] 

0 


(3 44) 


(3 45) 


where = 1-H 

Numerical procedure to determine expansion coefficients 
IS discussed in a later chapter. 


[3.12J (^i:.NEil/iL HEAT BALMOh 


A heat halance applied to the heat exchanger at any 
point along its length yields following relations for bulk 
temperatures* 

\av ~ ^ ^^2av “■ ^2av (3 46) 

Hence at z a 1 

H [52av (i) - ®say <3 ^7) 
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Expressing bulk temperatures in series solutions, the above 
equations can be manipulated using orthogonality condition 
Eq (A 35') and Eq (A 34) to yield following result 

= (1-H) Oo H (S) (3 48 

At the ends of heat exchanger the above equation becomes, 

T, (i/) = (1-H) 0^ + H 

lav ^ o 

^2av 

[3 13] HEAT EXOHAE'GER BEEICIENGY 

The efficiency of counter flow heat exchanger is 
defined as the ratio of actual heat transferred to that 
which would be transferred if the heat exchanger were 


(3 49) 
(3 50) 


infinitely long 

Actual heat transferred 


I^^IO " ^lav 


(3 51) 
(3 52) 


Heat transferred for infinite length 


== °i^i f^^lO - ^21^ 


(3 53; 


Ratio of preceding quantities yields 


H < 1 


e = 


H 


1 -■ T 


2av(0) 


(3 54) 


H > 1 


^ = ^lav(L) ^l^^”^2av(0)^ 


{3 55) 
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or using Bq ( 3 49 ) and (3 50) 


H < 1 


P 1 - H p , 

e = + 1 


(3 56 ) 


H > 1 


e = (1-H) Co + H 


(3 57) 


physically as length of heat exchanger increases (i e 
approaches infinity) efficiency must approach to unity 

A 

4Ience as L ® , e -»■ 1 Therefore from Eq (3 54) to (3 57) 


A A 


H < 1 as I -»■<» , H » 




(3 58) 


H > 1 as L 1 » 

T2a^(0)^ H-l/B , -v 1 


(3 59) 


The senes expansions for the end bulk temperatures may be 
written as, 


\av 


= °o + 


2 

-X./ L 




+ h (Bi.n (h)> ^ ] t3 60) 


+ «2,n (^2>’i 


(3 61 ) 
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Letting L -►» and applying Eq (3 58) to (3 6l) land 
Sq (3 59) to Bq (3 60), 

00 

H < 1 as 1 -> « , ^^2,n ^^2^^ 0 (3 62) 

H>1 asL-^“, -^n \ ^^i,n ^^1^^ ^(5 63) 

These equations reveal little about individual behaviour of 
the coefficients and 0^ as L However, computations 

have revealed that the terms in each of the above summations 
are identical in sign for all values of nj therefore, since 
summations approach zeio as L * 


H < 1 

Cn 0 

A, 

as L 

(3 64) 

H > 1 

A^ 0 

A 

as L « 

(3 65) 


[3 14] LONG HEAT BXGHAITGBR 


For long heat exchangers i e,, L " Eq (3*64) and 
(3*65) reduces series expansions for the temperature distri- 
butions of the heat exchanger fluids 


H < 1 


'^ 1 ^ ^1 * 


z) 






O ^ A 


(3 66 ) 
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. X * » . -|3 ^(L-z) 

(3,67) 
(3 67) 

-X 2 * 

H >1 i) + I ^n\,n ® 

n=± 


T2(^2. “ Oo + J, °n ® 2 ,n (^2’ ® 


(3 68 ) 


X 2 t 

•An 2 


n=l 


(3 69) 

Por fully developed region, only first term of the above 
expansions are required as other terms are negligible in 
magnitude compared to first term 

[3 15 ] fUllY LnVLLOPliD RBGIOI 

For fully developed region temperature distribution 
IS determined accurately enough retaining only first term 
in series expansion 


H < 1 ^pr^, S) ” e 


-a^®(£-z) 


(3 70 ) 


f Q ■”'1 2 1 '' 2 ' 


H > 1 T3_(nj_, z) =* Cq + ® 


2 * 

-^1 ^ 


(3 71) 


(3*72) 


2 ^ 

- . -h ^ 


^ 2 ^^ 2 ' 2) “ ^ 1 ^ 2,1 ^^ 2 ^ ® 
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Liquid metals have large thermal entiy lengths compared 
to nonmetallic fluids Pully developed region, for liquid 
metals may not te realized for small or moderately long 
heat exchangers Even for long heat exchangers, for accu- 
rate description of temperature at entry and exit, more 
number of terms are required 


CHAPTER 4 


PijUG PEO¥ lEDALIEATIOI 

[4 1] INTRODUCTION 

In principle assuming that sufficiently accurate 
values of f (r ) and ) can be found from recent expert— 

mental results for liquid metals, the Lq. (3 34) to (3 37) 
can be integrated with stated bound a^’y conditions Thus 
auxiliary functions P, P'*, G and G* may be obtained This 
may not be straight forward as the differential equations 
may not submit to analytical solution. Expansion coefficients 
will have to be integrated nume rib ally Obviously this would 
not desirable if it could be avoided* Much of engineering 
research can be associated with attempts to discover simpli- 
fications and convenient computational procedures 

Por turbulent flow of fluids, considerable simplifi- 
cations are offered if plug flow model is assumed Then it 
does not become necessary to resort to numerical techniques 
for determining auxiliary functions, in our problem 

[4,2] J’I'UG PLOW MODEL 

The plug flow model involves the assumption of unifbrm 
velocity profile i e , g^(r^) 1 and if in addition we neglect 

effect pf eddy diffusivity i e , fi(\^ ~ it is possible tp 



RELATIVE THERMAL RESISTANCE 


FIG 4 1 



RELAri/E IIISTANCE FROM WALL 


Relative t mal resistance 

' )lds 
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obtain analytical solution of the Bq (3 34) to Eq (3 37) 

With these approximations the equations for energy balance 
take Simplified form for region 1 and 2 In fact with plug 
flow and f 3 _(? 3 _) = 1, Eq (3 34) to Sq (3 37) take form of 
familiar Bessel equations and a relatively simple expression 
IS obtained for F, B , G and G^ 

Heat transfer analysis made with plug flow model for 
nonmetallic fluids give inaccurate results Stein (1966) has 
shown that with liquid metals this modelling yields meaningful 
and fairly accurate results, for low Peclet number (Pe ^ 100) 
Peclet number of 100 normally corresponds to s Reynolds number 
of 30,000 This range of Reynolds number for heat exchangers is 
definitely of practical interest Plug flow idealization is 
appropriate for liquid metals because steep temperature 
gradients are not localized near wall so that the error 
resulting from using uniform velocity profile 1 e , ^^(r^) = 
in this region is not significant For low Peclet number 
(Pe ^ 100) the assumption of f^(?i) = 1 is appropriate because 
eddy conduction plays a minor role compared to molecular 
conduction in the transfer of heat This is demonstrated 
by Fig ( 41 ) In our case Pe > 100 and hence an extension 
of this model will be needed This is discussed in Sec- 
tion [4 11] 
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[4 3] PLUG FLOW SCLUTIOUS 

With = ^ 1 (^ 3 .) = P» 37) hecome 


(f 


) + F (^1, X„) = 0 (4 1) 


ap(o. X ) 

= 0 , p ( 0 , x^) = 1 


A 

dr. 


4 - (?! 

dri 


dF*(r, , p„) 


dr. 


— > - h = ° 2 ) 


dF*( 0 , Pjj) 


dr. 


= 0 , P* ( 0 , = 1 


((rp + o) 

d ?2 


dG(rp, X„) 

I fa ' » ^ 

dr„ 


- 4 ? 2 +° ) ®(^ 2 > > 11 ’ = ° 


dG (1, X^) 

' ' = 0 , G (1, X^) = 1 


A 

drr 


, dG^Crp, p„) 2 2 


dr. 


( (r, + a) 


dr. 


) + w (rg+a) U (rp, = ^ 

(4 4) 


dG*(l, pj,) 


dr. 


0 t 


gMI, Pn^ 
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Solutions of Eq (4 1) to (4 4) in terms of P, F'* , G 
and G* are, 

F (r^, xj = Jo ^ ^n 5) 

P*(r3^, |3j^) = Iq ^ ^n 

^^^2’ ^n^ ~ ) (KjCo) X^(H- ff ) 

+ Ii[«Xn^l+a)] K:QCwX^(r 2 + cr)]) (4 7) 

(3^) = -f (1+cT) (J^[a)^^(l+a)]Y^[(.^^(r 2 +a)J 

+ Y^^Lco^j^d+o)] 

vhere , 

J^i Jq, Y^, Yq Bessel functions of first and 

second kind 

^1' ^ 0 * ^0 Modified Bessel functions of first 

and second kind 

Substitution of above expressions into Eq (3 30) and (3.31) 
yields eigenvalues Eigen functions can be obtained from 
Bq. (3 37) to (3 40) Finally expansion coefficients from 
(3 43) and (3 44) 
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[4 4j MEROW ARmjLUS APPROXIMillON 

A useful simplification is obtained letting R -»■ 1 This 
physically means a narrov annulus. As R -*■ 1 , a » We 
find that with R = 1, R and F"* remain unaifected Bq (4*7) 
an<i (4*8) can be simplified letting o and retaining only 
first term in large arguments approximations for 'Vessel func- 
tions A simpler approach is letting o » in Eq (4 3) 
and (4 4) which leads to following equations 


d^G(r2, X^) 


dr. 


S{J2. 


0 


(4.9) 


afj 



With initial conditions remainii^ unchanged 


0 


(4 10) 


It was found tl^at by both techniq.uss the expressions 
for functions G and O'* obtained were exactly seme Numerical 
computations have revesXed that for R as low as 0 5, error 
in computing efficieppy is less than 3.5 percent For 
present case R = 0 6^6, therefore we adopt narrow annulus 
app 3 X)Ximation approach The simplifications offered by 
employing narro'W annulus approach can definitely be used to 
sblve the problem, because the errors contained in the final 
results are insignificant 
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[4 5] auxiliary mug no ns 

For narrow annulus approximation the auxiliary 
functions obtained are as follows 


p 

x^) 

= Jo (’‘n 

(4 

11) 


^n) 

= lo (Pn> 

(4 

12) 

® (?2' 

Xn) 

= Gosh [ujXj^ (1 - r 2 )] 

(4 

13) 

G'Crj, 


= Oos [^^^(l-'rg)] 

(4 

14) 


For narrow annulus approximation the auxiliary functions are 
considerably simplified. Also R is eliminated except that 
it occurs in nondimensional parameter K 

[4 6] EIGENVALUE EQUATIONS 

Substituting F, F*^, G and G* in eigenvalue (3 50) 

and (3 31 ) we obtaih, 

«( (X„) = «> ‘“‘n 

sinh < «Di„) + K ooah (ol^) (4 14) 

and 

II (Pn) Bin (uPjj) + .Pn ^o^'n^ 
sin (.^) - KPr> ^'*•^5) 
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[4 7] N0EimiZIN& PAG TORS 

From equations (3 41) and (5 42) normalizing factors 
are given as, 


N. 


n 


sinh^ ((.x^) 


T o o sinh (2wX ) 

-i liK^ x/ Ji(x^) ] 


n 


(4 16) 


^n 


sxr? (|3j^) - (p^)] 


sin ( 2 w 6 ) 

5 Pa ii ( V ^ ] 


1 TTTr2 n 2 


n 


(4 17) 


[4 8] NORMALIZLD EIGBNFimOTIONS 


®l,n = 

- a)Xj^ sinh U\^) Jq ( X^^ 

(4.18) 

= 

“ Pn ^“Pn^ ^0 <Pn ?i)/V'IK'1 

{4a9) 

®2,'n = 

- KXn [u Xn^l“r2)]/V I\1 

(4,20) 

®2,n ~ 


(4 21) 


[4,9] EXPANSION GOBFFICIENTS 

In Eq (3 43) and (3 44) various factors for plug flow 


are given as follows 
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\ = 2/ K 


.2 o 2 __„2 / „ \ r 2 / 


Pn -Il (Pn)J 


Q ^ 




2* 

+ 6 1 


In."* 


n' 




1-e 

for n = k (4 22) 


2 sin(.«p^)][ ^ sin(o.|^)J[p^I^(P|^)I^(pJ 

- (^nh^V l^l]] 


for n k ( 4 23 ) 


%,k<^l)> = 2 / (J3_) ?3_ 


= - 2 [oip^ Sin (tt»p^)] [ sinh (“^i 5 .)][Pn ^n^ *^0^ 

+ \ lo (Pn) l»kl ) (Pn^ '4 24) 

®1 *®l,k ®l,n ^ ^ i ^l>n '^’^1 


- 2 [u)|3jj. Sin (u^)] [wX^ Binh (uX^)] [^ 

+ l»kl (4.25) 
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B 


1 \,k: = 2 / (f^) (^l)^■l dr 




+ e 


2^^ 

-Kh , .. ^ 

/(l- e ^ ) 


-A„^L 


for n = k 


(4 26) 


= 2[u,x^ sinh (u>x^)] Uxj^ sinh ((aX^^)] 


- Jo( V3 /[(VIhJ IhJ )(Xj^2 _ 


for n k 


(4 27) 




= 2 / (f,) 


= 2 




(4 28) 


(y|<i ) p„ 


®1 <\,n (%» = 2 / ®l,n (? 1 > ^1 '3^1 


oi\ sxnh (u)X„) J, (X«) 


= - 2 


n 


S — ■ I . . .C. S ... (4 29) 


(ir|if„i ) V 


RELA'^tVE THERMAL PES’STANCE 


FIG 4 2 



mal resistance distribution 
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[4 10] BULK TEMPER ATURE 

Expressions for bulk temperatures are given as, 


’^lav (2) 


wsin ( 0)6 )I ( B ) -e^(L-z) 


= G + I [2 A. 




n 


n=l 


n 


flH'l 


- 2 C 


o)Sinh iui\) 


n 


f|N, 


n g n ] (4 50) 


n 


’'2av (-) = 


°o + I [■V 

n=l 


sin(o)j3^) 


0) f|R^ 


K sinh (o) X„) J. ( X^) -X 




n 


(0 


f |N. 


] (4 31) 


n' 


It can be verified that the orthogonality condition is 
satisfied for plug flow, narrow annulus (R = 1) case General 
heat balance is also satisfied 

[4 11] PLUG PLOW EXTENSION EOR LARGE PBCLBT NUMBER RANGE 

Li(luid motal heat exchanger applications to nuclear 
tecnology involves very hi^ flow rates Conseouently 
Reynolds numbers are large (Re No ^10^) Peclet number 
corresponding to such range is much greater than 100, 

Por this range of |eclet number, in plug flow idealiza- 
tion, assuming = 1 is inaccurate For Reynolds number 

of 10^ there is a large variation in eddy diffusivity across 
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duct crosa-section as is seen in Fig (4 2 ) stem (1966) 
suggested an approximation to extend the validity of plug 
flow for peclet number range up to 1000 This approximation 
involves the use of a properly averaged space independent 
total conductivity which results when a simple variable 
change is made 


The following derivations follows the method suggested 
by Stem (1966) 

[4 12] BhHIVATION 


A variable change in radial variable r^ is introduced 
The new variable is related to r^ by a first order 

differential equation 


Tube 




J ^ V 

fl(rp 


(4 52) 


Annulus 



+ cr 
r2 + 0^ 


+ 





(4 33) 


•• 

with the boundary condition x^( 0 ) = 0 and with a constant, 
so oho sen that x^(l) == 1 
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Introducing variable change in the equations (3 12) 

and (3 13)} in terms of they can now be written as 
follows 


lube 


1 3 







, 2 ) 
dZ 


(4 34) 


Annulus 


(x2 + cr) 3X2 


[(Xg + a) — ^ ] 


3 X 2 




•Un ^ 


ai^Cx^jz) 

i 

m 

3X, 


whe re , 


z S5 kj z j and L =s I 


A , A 

r, p drT 

= g;]^(r3_) (j-) — 

^1 °^l 


+ 0^2 <3r2 


gpixp) = S„(r2) (~ ) 

^ ^ X, + a 


dx. 


(4.35) 


(4 36) 
(4 37) 


(4 38) 


The boundary conditions become, 
Entrance ., 


(x^, 0) 


0 


(4.39) 
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iM M 

^2 (^2, L) = 1 

Interior 

3 ( 0 , z ) 

^ 0 


(4 40) 


t4 41) 


K 



3^1 (1, z) 



312 2 ) 

3 X 2 


sIt ( 1 , z) 

K;^ ^ T (0,z) - T (l,z) 

3Xj^ 


sL (l,z) 

— i = 0 

3X2 


(4 42 ) 


(4 43) 


(4 44) 


For turbulent flow and all ranges of Pe No the gCx^) and 
g(x 2 ) as defined hy Bq (4 37) and (4 38) have average value 
over cross-section equal to unity For small Pe No g(x^) 
and £(^ 2 ) are approximately equal to unity for a large portion 
of cross-section of the heat exchanger duct 


liquid metals are less sensitive to the exact represen- 
tation of velocity profile compared to eddy diffusivity 
representation, Stein (1966) Therefore as an approximation 
it will be assumed that g(x3_) and g(x2) are identically equal 
to unity for all values of Peclet No. This reduces Bq (4 34) 
and (4 35) to following 
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Tutie 


t 





z) 


Annulus 


(^i, z) 
dz 


(4 45) 


r-i- V[(22 + 

+ a 3^2 


3T (x , z) 

a) — M 

3X2 



3T2(x2> 2 ) 


3X 


2 

(4 46) 


with boundary conditions remaining unchanged 

Bq (4 45) and (4.46) are lusntieal to Bq (3 12) and 
(3 15) when plug flow model 1 e , « 1 is applied 

to later set of equations, and maicing following changes in 
non-dimensional parameters 


kt 

K — K 

(4 47) 


(4 48) 

z ^ z 

(4 49) 

In effect, the k'^ approximation has changed molecular oonduc- 

«L> 

tivity to ill equations Hence the results 

obtained for plug flow without k^ are valid for large Pe Ho 


approximation with above modification in parameters only 
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[413] Determination of Ic^ 


eqns 


Expressions for may "be determined from integrating 
(4 32) and (4 33) between 0 to 1 




1 

f 

0 

dx^(^l) 
il (rj^) 

1 

f 

1 % 

1 

^1 ( ^1 ) 

1 

dX2 (r2) 

i 

(X2(r2)+ O') 

1 

( 

1 

/ 

0 

1 . 



(4 50) 


(4 51) 


There are two complications m determining k"^ First it is 

necessary to know fi(r^) which essentially is knowledge of 

£ as a function of radial location This can be obtained 
Hi 

from experimental results Second complication is that right 
hand sides of equations (4 50) and (4 51 ) take indeterminate 
form at r^ =: 0, i e , x^ = 0 

An alternative method is suggested by Stem (1966), 
according to which liquid metal heat transfer coefficients 
can he computed from plug flow values by simple relation 


(b) = (2)3pj' 


( 4 - 52 ) 
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where, 

PP - Plug Plow 
Nu - Nusselt Number 

Since detailed temperature distributions are not reqiuired 
there is no need to determine f^rom equations (4 32) 

and (4*33) Appropriate values of can be determined 
from equation (4 52) when applied to fully developed 
nusselt numbers either from experiment or by computation 

The equations used for evaluating k^'*' are , 

Tube 

k^"^ s (7 + 0 025 Pe° ®)/8 0 (4 55) 

4nnulus 

kp'^ =» (5 25 + 0 D188 Pe^ ® (^)^ ^V6 0 (4 54) 

^21 

Sq (4 53) is obtained from oorreljatio ns given by Lyon (1955) 
and Bq ( 4 . 34 ) from correlations given by Bailey (1950) 


CHAPTER 5 


RESULTS AHD DISGUSSIOH 


[5 1] El GEE values AND EIGBNPUUCTIONS 

P 2 

The square roots of the eigenvalues, and , 
were determined from the roots of the Eq (4 14) and (4 15) 
These nonlinear algebraic equations involve nondimensional 
parameters R, and The values of these and other 

parameters for the Kalpakkam heat exchanger are given helow 
To calculate these values thermal conductivities of sodium 
and austenitic stainless-steel were obtained from Raznjewi 
(1976), 


Pe^ 

S3 

160 


js: 

190 

V 

=s 

1 055 

2 

ss 

1.11 

K 

as 

( 5722) 

H 

9S 

1.65 


)s 

( 5698) (kj_+) 

#v 

z 

s 

(2 6)(k;i^+) 

The values 

of 

kj^ and obtained 


and (4.54)* Numerical values given within 


from Bq (4 55)» 
the brackets are 


plug flow values 
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llie Bq ( 4 14 ) and (415) having J<, and 00 as 
multiplying constants were solved by Regula-Palse method 
The accunacy in determining the roots of these equations 
IS of the order of 10 Once the eigenvalues are known 
it is easy bo calculate eigenfunctions f 2 X)m Eq (4 18) to 
(4 21) 

[5 2] expansion GOEPPIOIBKTfa 

The expansion coefficients were calculated from the 
Bq (4 22) to (4 27) and the nonhomogeneous part from 
Bq. {4 28) and (4 29) To obtain and an infinite set 
of equations must be solved It has been shown by Brown (1968) 
that it IS enough to consider only ten equations to obtain 
Oq accurate up to three decimal places 

Bor long heat exchangers i e , z ^ 1, only first term 
was required to calculate temperature distribution The subse- 
quent terms had negligible contribution except close to the 
entrance or exit of the heat exchanger 
[5 3] COMPARISON OB PRESENT MBTHOP WITH NTH-METHOB 

The performance of a conventional heat exchanger is 
usually predicted by the effectiveness - ITU method In this 
method, the calculation of NTU (Number of heat transfer 
units) involves the estimation of heat transfer coefficient 
for primary and secondary side, using 

TT^rsfTTi 

Acc* No vs «.*<»!» 
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The UH]? approach over predicts efficiency, but the 
error is not serious because the actual efficiency quoted 
by manufacturer is close to hundred percent It is also not: 
very clear that manufacturer's data represents experimental 
values or calculated values Table 5 1 also sho-ws that UHF 
over predicts exit temperature for secondary sodium and 
under predicts exit temperature for primary sodium The present 
method shows the opposite trend Temperature variation for 
UHT' and present method along length is shown in Pig (5 1) and 
Pig (5 2), The maximum temperature difference in two 
methods are at the outlet since the inlet temperatures were 
assumed to be same in both cases 

Manufacturer's quoted temperature values are only 
available at present Actual temperature values will be 
known only after experiments conducted at Kalpakkam These 
experiments might throw light upon the actual obtainable exit 
temperature values At present the Inteimediate heat exchanger 
IS under fabrication by BHEL 
[5 4J iiPPhOT OP VABIATIOE OF PLOW RATE 

Under noimal operating conditions for IHX, H si 65* 

A change in operating reactor power level wiH have a corres- 
ponding change in the quantity of heat transferred at IHX 
from primary to secondary sodium. If reactor power 

level goes up more heat will be exchanged at IHX Consequently, 
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a higher secondary mass flow rate will be required to 
transport this heat If reactor power level goes down 
less heat is exchanged and a lower secondary mass flow 
rate will be required In either case H will vary as it 
represents ratio of primary sodium to secondary sodium 
heat capacity flow rates 

It IS reasonable to assume that primary flow rate and 
its exit temperature from the core is constant ¥hen power 
level goes up, the higher power is transferred to primary 
sodium inside core Ihis can be achieved keeping same 
mass flow rate and exit temperature, from reactor core, with 
a lower inlet bemperature to the core When power level 
reduces, other conditions remaining constant, a higher inlet 
temperature to the core is anticipated 

Secondary sodium mass flow rate will depend upon the 
quantity of heat to be transferred But it is reasonable 
to consider exit temperature of secondary sodium from steam 
generator to be constant This is possible because a 
changed secondary mass flow rate will result in a change in 
the quantity of steam generated for same exit temperature 
from steam generator Only IHX exit temperature will 
increase or decrease depending upon increase or decrease in 
reactor power level This may be of practical interest also 
if in secondary loop due to blockage or additional resistances, 


FIG 53 
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the secondary sodium mass flow rate changes Any change in 
secondary sodium flow rate due to either reason, will result 
in a corresponding change in H With the reasonable assump- 
tions of keeping inlet temperatures to the IHX constant, 
exit temperatures for both primary and secondary sodium 
can be predicted, under steady state conditions S'lg (5 5) 
and Pig (5 4) gives such a plot, both for present method 
and UHP method 

As anticipated for higher values of H both predict 
hi^er exit temperatures for secondary sodium Por H decreas- 
ing both methods predict decrease in secondary sodium exit 
temperature This trend is observed for primary sodium also 

It IS interesting to see that for H > 1 UHP over predic 
the exit temperature values for secondary sodium and under 
predicts for primary sodium, which essentially means over 
prediction of efficiency compared to the present method Por 
H < 1, It under predicts efficiency A good agreement 
is obtained for H s= 1 This means that for H = 1, UHP method 
and present method yield same results 

Por H> 1 and increasing H will decrease which 

for H > 1 is given by secondary heat capacity flow rate 
Sxnce NTU is defined as decreasing will 

increase HTU (See Appendix B) and from Eq (B 6) it is obvioi 
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that efficiency will increase This increase in efficiency 
IS much higher for UHP compared to present method for equal 
increase in H, Hence UHP method over predicts temperatures 

Por H < 1 and decreasing H will fix U ^ hecause now 

min 

^min ‘^o^^sponds to primary heat capacity flow rate which 
IS considered constant Prom Eq. (B 5)» it is clear that 
denominator does not decrease as fast as numerator Conse- 
quently efficiency increase for H < 1 and decreasing H is 
at a slower rate than for H > 1 and increasing H This 
resulted m under-prediction of efficiency for H < 1 by UHP 
method compared to present method 

Por H w 1, for a heat exchanger energy balance equa- 
tion at any point along the length predicts constant bulk 
temperature difference between primary and secondary 
sodium This means that there is a constant heat flux all 
along the length Therefore UHP should predict accurate 
results under such condition Present method also yields 
same results and therefore it is advisable to use NTU method 
when Hal, based on UHP conditions 

[5 5j COMPARISON OP WITH PLUG PLOW MOBBL 

The Pig, (5 1) to Pig (5»4) were presented using k^'** 
approximation. Neglecting k^"*" in the present case results 
in the under prediction of primary sodium exit temperature by 
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3 0 and over prediction of secondary exit temperature sodium 
by 4 0 The Fig (5 5) shows predicted temperature profiles 
along heat exchanger length for the two cases, viz , with 
and without approximation The peclet numbers for the 
present heat exchanger are not very for removed from 100 
Therefore, the effect of eddy conduction incorporated using 
IS not very pronounced The two temperature profiles 
are therefore very close to each other except at the ends 
where marimum difference is observed 

[5 6] VimiATIOh OF HEAT FiUX ALOm HBAl LXCHAiNTShR LEMIH, 

The heat flux for the present heat exchanger varies 
strongly along length From Fig, (5.6) we see that the 
heat flux at one end is ten times ab the other end Hence it 
IS incorrect to assume uniform heat flux boundary condition 
for present heat exchanger But results obtained for the prese 
case even assuming UHF do not lead to serious errors This is 
because for H > 1 the over-predicted value of efficiency is 
near unity and quoted value is also very close to unity 

For short heat exchangers with smaller efficiency and 
H 9^1, UHF assumption may lead to gross over prediction or 
under prediction of the heat exchanger performance But when 
H « 1, the present method may yield results which may be very 
close to UHF predictions 
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[5 7] liAHAb TEi''IPiiEATUEE VARIATION 

UHP meuhod can predict only bulk temperature variatio 
of th.e heat exchanger along the lengch But present method 
was evolved solving fundamental eq^uations for energy balance 
The solution thus gives temperature at all points within the 
considered boundaries Hence present method provides more 
information bhan TJHB method 

A detailed knowledge of temperature at all points, 

IS of interest to the designer The radial temperature 
profile gives the magnitude and location of hot spots The 
Fig (5 7) shows that at the centre of heat exchanger, 
maximum temperature occurs at the outer boundary of the 
annulus The present method can be used to determine such 
radial temperature profiles, at any distances from either 
ends* 

Only first term is req[uired to calculate the radial 
temperature profile in fully developed region, i e , in the 
middle of heat exchanger,. But close to the entrance and 
exit more number of terms are required Number of terms 
depend upon the accuracy with which temperature is to be 
predicted Foi short heat exchangers fully developed region 
may not be realized at all This requires large number 
of terms to calculate radial variation along the length 
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Prom Pig (5 7), it is found that the temperature drop 
in tube metal is maximum, i e , 17°C Temperature drops in 
primary sodium is ll'^G and in secondary sodium it is 16 5^0 
Temperature drop is proportional to the resistance offered 
by the medium Prom Pig (5 7), we conclude that for the 
present caso resistances in all three regions are of the 
same order of magnitude 

[5.8] Sa&GL&TLh fUTUaeJ WOBK 

a) Recent eddy diffusivity and velocity distribution 
correlations for liquid metals can be made use of in deter- 
mining accurately It was not required for the present 
problem because the Peclect numbers were not much greater 
than 100 

b) The present method involves use of plug flow model 

where we assumed f]_(? 3 _) = ^ This assumption 

can be removed if experimental correlations available for 
flow along rod bundles and flow inside tube are employed to 
determine and The basic energy equation 

can then be solved numerically as a two point boundary value 
problem Such a general solution will be of immense value 
It can be also used to ascertain validity of experimental 
correlations and validity of any approximate theoretical 
approach like equivalent annulus approach and RTF method 
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APPEI^DIa a 


[A 1] SEPAEATIOiM OP VAPIABLES 

The solution of Bq (3 12) and (3 13) involves H, K, 
Ky and L as nondimensional parameters The solution 
requires the knowledge of 

determined both by Prandtl and duct Reynolds number Values 
of these nondimensional parameters determine for a specific 
heat exchanger, what may be called as operating conditions 
for gheat exchanger 

To solve Bq (3 12) and (3 13) separation of variable 
method is adopted Let, 

Ti(r^, z) = ®(2) 1) 

substituting Bq (A 1) in the Eq {3»12) and (3*13) results 
in Eq (A 2) and Eq (A3) 


Tube 


1 

^ . A. 

e( z) 


d6(z) 

A 

dz 





drn 


(A 2) 


Annulus 



de(z) 

ii “ g2(?2) Eg ir^)(}2+a) dJj 


^ dBpCr-) ^ 

[(^2 + d ) ] 


dr. 


(A 3) 



ts> 


A- 2 


The left hand sides of the above equations are functions of 

A 

only whereas the right hard sides are functions of 

only , consequently, they can he equal only if they each 

2 

equal a constant, - X 


de(z) 


- X' 


(A 4) 


e(S) d(z) 


2 ^ 

\ -X z 

or e(z) = e 

The right hand side of Eq (A 2) and (A 5) can he 


(A 5) 


dr. 




(r, ) 2 / * V * -cl N c 


dr. 


( A,6) 


and 


^ dE (rp) 

— [f2(r2)(r2+a) ;; ] 

drp 2 ^ ^ drp 


- ru (r2+a)E2(r2) 

(A 7) 


The intenor boundary conditions become, 


dh^(O) 


= 0 


dr. 


K 


dE^(l) 


dr. 


dSplO) 


dr. 




dB^(l) 

d^“ 


B^iO) - B^(l) 


(A,8) 

(A.9) 
(A 10) 



0 


A-3 


A 





(A 11) 


[A 2] TWO BA&IOli bTUEiyi-IIOUVILlA PROBW''! 

The Bq (/ 6) and (A 7) correspond to -what may be 
called a two region Sturm-Liouville problem and may be 
written as a single differential equation 

(k ^) + ^2 gy = 0 , 0 2 (A 12) 

with boundary conditions, 

y' (0) = 0 and y' (2) = 0 (/,13) 

where x = 0 and x = 2 correspond to r^ = 0 and 1*2 = 1 
respectively 

The functions y, k, g are defined over both regions, 
tube and annulus, of the heat exchanger by. 

Tube 


(x 

A 

^ ^1 

\ 0 < 1 1 ) 


y 

= E3_( 

:^i) 

(A 14) 

k 

A 

= ^1 


(A 15) 

g 

A 

= ^1 


( A 16 ) 
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Annulus ( v 

^2 ^ » 

0 1 ?2 — 


y == 

E2(r2) 


(A 17) 

K 

(r2 + o) 


(A 18) 

g 

- to (r 2 

+ a) ^ 2 ^X 2 ) 

(A 19) 


where the function y must meet tne compatihility conditions 
at the inner wall of heat exchanger Eq ( A 9 ) and (A 1^) 
and the functions k and g have finite discontinuities at 
the inner wall of heat exchanger 

The difficulty in applying the £>turm-Iiouville theory 
to the present system lies in the fact that it specifies 
that the functions must be continuous The functions are 
not continuous in the system being considered Because of 
the fact that the functions are continuous in their regions, 
1 e , piece-wise continuous, Stein (1966) has been able 
to show the Sturm-Liouville theory applies to the present 
system 

According to Stein, then, there are infinite sets of 
eigenvalues and corresponding eigenfunctions, having limit 
points from - « to ® Thus in one region problem the 
negative eigenvalues do not exist but for two region 
problem, having discontinuity at interface, both negative 
and positive eigenvalues exist 



The Eq (A 4), (A 6) and (A,7) according to above two 
region bturm-Liouville problem have infinite sets of positive 
and negative eigenvalues and corresponding eigenfunctions 
These equations thus may be written as, 


z) _ .2 

e(a) dz - ^ 

^ OT 6 CZ-) =. £ 

d / A V A y, ( ) 


(A 20) 
CA Zi) 




(A 22) 


4- [fj (r*2) 

dr2 ^ drg 


- (Jg+ct) g2(J'2) B2,n^^2> “ ^ 25 


and boundary conditions; 


dE (0) 


(A 24) 




(A 25) 






(A 26) 


^^2.n 


(A. 27) 


, +«B 
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A-6 


aAr 

where is the eigenfunction corresponding to 

[A 5] SOLUTION 

The general solution to the Bq (A 6 ) and (A 7 ) may 
he written for region 1 and 2 as follows 


Tube 


A A. 


n = + » 


1 2 * 




(A 28) 


Annulus 




Iir. 


2 * 


(A 29) 


where 0^ a Expansion coefficient associated with n-th 
eigenvalue 


In order to determine eigenvalues, eigenfunctions and 
expansion coefficients we first of all derive orthogonality 
property 

[A 4] ORTHOGONALITY OP EIGBNPUNCTIONS 

The orthogonality property of eigenfunctions may be 
derived in the following way Bq (A*22) is first multiplied 

A 

by B, „ to obtain a relation between different order 

J* ^ Ul 

eigenfunctions A second relation is obtained by inter- 
changing m and n in first relation. The two expressions 
are then integrated between 0 and 1, subtracted and 
simplified to yield Bq (A 50) below 


A-? 


(X 


n - (ri)> 


- 2E. 


. n, (^) * v,(l) 

( 1 ) [ Bi.eCi) — ]( a 5 o) 

drn dr^ 


A similar treatment to Eq (A 23) yields Bq {A 31) given 
below 



2 A dE2 Y,(0) dE- (O) 

= IT [1^2, m^O) ^ 51) 


where n = + 0,il,+2, 

The Eq (A* 30) and (A 31) are combined by applying the 
innerwall boundary conditions to yield 


A A 


( Xn " Xjn ) [B3_ t%^n(ri) 


- HB2 ^ 


(A 32) 


m 

~ i 0, + 1, + 2, 

n ~ 

This gives desired orthogonality condition 


Bi - HBj (Jj) E2,„(r2)> = 0, 

n m (A 33) 


A-8 


- HB, {[B. (J ]2) = 


2n'^2-J ^ \ f n = m 


(^ 34) 


n 


i 0 ) i 1 > + 2 , 

where { } and B2 { ) physically mean average over cross- 
section They are defined as follows 

Tube 


Bi { } = 2 / gj^ir^){ } r^ dr^ 


(A 35) 


Annulus 


Bg { > 


1 + 2o o 


•j* 

I 82(^2^ { } (rp + a) d^p (A.36 ) 


Thus the average temperature over the cross-section (bulk 
average) is given as { T^ ^ » which is denoted as 

\av ^ 

Note that dimensional equivalent for bulk averages 
are defined for tube and annulus, as follows 
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;i2 12 

s i2n f ^ }dr^)/(27c dr,) 

0 0 


(A 37) 


A-9 


Annxilus 

^22 ^22 

^2^'^ = ( 271 J U2(r2)(r2){ } dT^)/{2Ti dr ) 

^21 ^21 ^ 

( A 58 ) 

[15] BIGLi\[OTGTIONS 


Solution of the Bq (l 22) and (/ 23) for = 0 

are S]_,o ®2,0 equal to a constant No generality 

IS lost in assuming those constant equal unity 


A A 



The zeroth order nomalizing factor as given by Eq 
then becomes, 


(A 59) 
(5 61 ) 
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= l-H 


( ' 40) 


solution of Eq (A 22) and )(i 23) is best effected by 
defining auxiliary functions ^j^) and ^(^2* ^n^ 
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where P{r 2 _> X^) and ^(^2* solutions to 


(A 41) 
(A 42) 


dr. 


[fl(q) ri 


dl’C r. 




dr. 


] + X 


n 


g3.(r 


l) P(rj_, 




(A 43) 
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and, 


^ [f2(^^) (V-’) — f. ■— J - 

*^^0 dr^ 


2 . 2 
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^ ^^2' ^n^ ( A 44 ) 


XQ ijiJ.aI-QQnditiQns 


di(0, X^) dF(l, x^) 

jj = Oj 

dr. dr„ 


= 0 


F(0, \^) = 1 j G(l, X^) = 1 , 


n' 


( I 45) 
(A 46) 


The constants are so defined that they satisfy inner vail 
boundary condition Eq (A 25), 
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dB. (1) 
l.n^ ^ 

A 

dr. 


^^2.n (°) 
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(A 25) 


1 “"-2 
which in terms of auxiliary functions become, 
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( iA 47 ) 
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Therefore let, 
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{A 48) 
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y) 

(A 49) 
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where 

n ss 

Hh 1, + 2, + 5» 

(A 50 ) 
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[A 6] J*IG-£NVALUL EQU/TIOH 


iho eigenvalues are determined by the innerwall 
boundary condition Eq (A 26) which in terms of auxiliary 
functions P and G- becomes 




as(0, di'd, x^) 


di , 
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dF(l, Xj 


n' 
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G(0, X„) 
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dG(0, X ) 

— P(l, X^) (A 51) 


A 
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^ ~ ± 1» ± 2, + 3, 


The eigenvalue equation is therefore defined as, 



dG(0, X ) dP(l, x„) (iG(0, X„) 

IS, — JT-S -:r-^ + — F(1. X„) 

dr^ dr^ drg 


The eigenvalues 
equation, 
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(/ 52) 


are obtained fxom the roots of the 


(♦ ) - 0 
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[A 7] iiXPMSIOI GOJSiFIGIENTS 

Temperature distributions in tube and annulus are 
given by Eq (A 28) and (A 29) 
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Relation defining ejjpansion coefficients, 0^, are obtained 
from tlie entrance and exit temperature profiles, of the heat 

exchanger fluids 
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The end temperature distrihutxons may be written m terms 
of the senes solution, 
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(A 58) 
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J '■2,n 


(A 59) 
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2 

- X I 


= I °nh,n " " 

n=: - ® 


(/ 60) 


■(- n> 


n=-« 2,n 2 


2 '*■ 
- X ^ L 
n 


(A 61) 


Manipulation, of above equations yields expressions for the 
expansion coefficients Bq (yi 58) and {a. 59) are multiplied 
by and respectively^ and the bulk average operation 

is applied to the result The expressions are then subtracted 
and the orthogonality condition, Bq (A 33) and (A 34) applied 
to give Bq. (A 62) below A similar treatment to Bq ( 60) 

and (A 61) yield Eq (A 65) 


\ = - HBj (rp P2(r2)> 


(A. 62) 


” T1 

0 W e ^ 
^n ^ n 


, 2 r 

X„ L 


A A 




-HBj (rj)! 


(A.63) 
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the above expressions cannot be used directly as they 

(rg) However, 
these unknown functions can be eliminated by substituting 
thorn as their series expansions After substitution the 
two Bq (a 62) and (A 63) are added to give Bq (A 64) 


contain the unknown expressions Pq(:^^3_) and Pg 





N. 


n 


(1 + 




^ ®1 ^^l,n ^^1^ %,k 


’ HB2 CE^^^tr^)) 


(A 64) 


n 


0 > ± + 2 , + 3 , 


For n 0, tho term for K == 0 drops out of above summation 
and Oq does not appear in the expression An explicit relation 
for in terms of higher-ordered expansion coefficients may 
be obtained for n = 0 


Foi n = 0, Eq (A.64) becomes, 


A 





k = -» 
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(e 


2*- 


®1 ^\,k ^^1>^ 


- H 

(A 65) 


A fairly straightforward manipulation of the Bq (A 64) 
yields the following simplified form, 
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ii 


6 , 

n 
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1 { l~e 
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(A 66 ) 


where, 


*n,k = “l,n \,k “ (^67) 
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n,n = 


2 '" 


1-e 




-) \ ] j k = n 


(/ 68) 


^ = t 1» t 2, + 3, 


The zeroth order expansion coefficient does not appear in 
the atove equation and must he calculated by Eq (A*65 ) 

[A 8] NOEMALIZEE EIGENFUHGTIONS 


At this point it becomes desirable from practical 
viewpoint to normalize the eigenfunction with respect to the 
square root of the absolute value of Absolute value is 

used because for n < 0, becomes negative This normaliza- 
tion is desirable for computational considerations, since 
it was found that unnormalized eigenfunctions become 
unmanageably large as their order is increased Let 



(A 69) 



^2.n ^^2^ 


/!\l 


(A 70) 
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'll! the sigenfunc tions are normalized except zeroth—order 
e igenfunc tions , which remain, 

® 1,0 ® 2 , 0 ^ 

using normalized eigenfunctions, the orthogonality condi- 
tions become 




B. 


n ^ m 
2 


1 (rj)]'^). 


(A 71) 
+1 n > 0 
-1 n < 0 


n 


(A 72) 

Relations which were derived so for one still valid with the 
value of taken as + 1 for n > 0 and - 1 for n < 0 

[A 9] BhjFIinTIORS TO AVOID USE OR NEGATIVE INDICES 


Eigenvalues are given as 




X 


2 

m-1' 


X 


2 

m 


""00 < m ^ ® 

2 

We define auxiliary eigenvalues 


where 


Such that 


A -n 


(/ 75) 


Eigenfunctions associated with negative indices are 
also redefined as, 


ii.. 


-r, = 


1,-n 


^l.n (-1) 


(/ 74) 


^2,1-n ^^ 2 ^ 


^2,n ^^2^ 


(A 75) 


Expansion coefficients associated with the negative eigenvalues 
may be redefined in terms of auxiliary expansion coefficients 


and 


0 -►0 
n n 

2 * 
-e I 

«.n - 


(A 76) 
{ 'I 77) 


2 2 

The above definitions now give summation for and 

from 0 to + * and thus obviates use of negative indices 


APPi^roiX B 


HEAT TLirX SiMx'LE OALGJLATION BOR 
HAlPAKi£AM DATA 


[B 1] NUMBER OB HEAT TRAi^SBBR IRilTS 

The number of heat tratisfer units (RTU) for a 
heat exchanger is defined as "where UA for double 

pipe IS defined from equation (B 1), 


1 

UA 


2'n;k..L 

w 




2k. 


"w 


)i2^2 


+ in (^)] 
^22 


(B 1) 


and C =! Minimum heat capacity flo"w rate 
min 

A^ = Surface area of tube 
A 2 = Inner surface area of annulus 
»: Conductivity of tube metal 
= Heat transfer coefficient for tube side 
= Heat transfer coefficient for annulus side 
= Length of heat exchanger 


^1 

h2 

L 


For double pipe heat exchanger tke following 
convenient expression for OTU, in terms of Husselt numters 
can be derived from Eq (B 1), 




min 

2nk^L 


2k 


w 


Nu^k^ 


{^22-^Pl^ 

Fu2 k^ 



(B 2) 


The Nusselt numbers were obtained usmg experimental 
correlations of Loyn (1955) for tube and Bailey (1950) for 
annulus 


Tube 

= 7 + 0 025 Te^° ^ (B 3 ) 

Annulus 

Nug =« 5 25 + 0 0188 ® (— )*^ ^ (B 4) 

^ ^22 

In Bq (B 2) various parameters for Ealpakkam heat exchanger 
have following values 


Bu^ a 44 

Nu2 « 6 676 

= 19 77 ¥/m°C (at 425°G) 

= ^2 watts/m°G )(at 425^0) 

L « 2 48 m 

^21 *= H X 10”5 m 

^12 = 7 X 10*“^ m 

and ^^^22*"^21^ ** Hydraulic diameter 

4 3c Gross -sectional area of Annulus 

) 


( 


Wetted Perimeter 


B 3 


With these yalues Eq (B 2) gives 

JTU = 10 470 
[ B 2] EPPIGIENGY 


Efficiency for a heat exchanger is defined as the 
ratio of actual heat transferred uo that which would he 
transferred if heat exchanger were infinitely long Por 

counter flow heat exchanger, efficiency interms of ITTU 

nnd TT Ch - Capacity Plow Rate . 

Secondary heat Capacity Plow Bate ^ given as 


H< 1 

1 - Exp (-NTU (1 ~ H)) 

e (B 5) 

1 - H Exp (-NTU (1-H) ) 

H > 1 

1 - Exp (- BTU (1 - -i)) 

s £ (B 6) 

1 - “ Exp (-WTU (1 - -^)) 

Prom equation (B 6) efficiency can be calculated suhsti- 
tuting the values of NTU and H Por present case H =1 65 
Therefore Bq (B 6) gives 


e 


0 993 


[B 3J OBTLj^O} TEMPBRATUliJib 


B h 


1 


I 

1 



Efficiency for a given heat exchanger in terms of 
temperatures is defined as follows 

H>1 


e =: 


^IL ~ ^10 
^21 - ^10 


(B 7) 


H < 1 


^ ^2L ^20 

e = T — 

21 no 


(B 8) 


where subscript 

1 denotes tube side and 2 denotes annulus side, 

0 denotes entrance side and 1 denotes exit side 

From general heat balance we can write, 

(hx, - ho) = H (Igi, - Tjo) (B 9) 


For present case secondary outlet temperature, can 

be determined from Bq (B 7) liimary outlet temperature, 
T 2 q» can be determined from Bq (B 9), They are given 
as , 
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